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Abstract. The Daya Bay experiment was designed to be the largest and the deepest underground among the many current- 
generation reactor antineutrino experiments. With functionally identical detectors deployed at multiple baselines, the experi- 
ment aims to achieve the most precise measurement of sin 2 2613. The antineutrino rates measured in the two near experimental 
halls are used to predict the rate at the far experimental hall (average distance of 1648 m from the reactors), assuming there 
is no neutrino oscillation. The ratio of the measured over the predicted far-hall antineutrino rate is then used to constrain the 
sin 2 2613. The relative systematic uncertainty on this ratio is expected to be 0.2~0.4%. In this talk, we present an improved 
measurement of the electron antineutrino disappearance at Daya Bay. With data of 139 days, the deficit of the antineutrino rate 
in the far experimental hall was measured to be 0.056 ± 0.007 (stat.) ± 0.003 (sys.). In the standard three-neutrino framework, 
the sin 2 20]3 was determined to be 0.089 ± 0.01 1 at the 68% confidence level in a rate-only analysis. 
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INTRODUCTION 

As fundamental particles in the standard model, 
(anti)neutrinos were initially thought to have zero 
mass. Such an assumption was supported by the exper- 
imental evidence that only left-handed neutrinos (also 
right-handed antineutrinos) were detected [1]. However, 
in the past decades, the phenomenon of neutrino flavor 
oscillation observed by Super-K, SNO, KamLAND, 
MINOS, and many other experiments successfully 
established the existence of non-zero neutrino masses 
and the neutrino mixing. A recent review can be found 
in Ref. [2]. The neutrino oscillations are commonly 
described by the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix and two neutrino mass-squared differ- 
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ences (Ara^ := —m^ and Arajj := m\ —mf) [3, 4, 5]. 
The PMNS matrix denotes the mixing between the 
neutrino flavor and mass eigenstates. It contains three 
mixing angles 612, 623, and 613, and an imaginary phase 
8, referred to as the CP phase in the leptonic sector. 

As of two year ago, 613 was still the least known 
among all three neutrino mixing angles in the PMNS ma- 
trix. The best constraint was from CHOOZ reactor an- 
tineutrino experiment with sin 2 20i3 < 0.17 at 90% con- 
fidence level (C.L.) [6, 7]. A global analysis [8] in 2008 
including both solar and reactor neutrino data suggested 
a non-zero 613. However, in the past 18 months or so, an 
explosion of data from multiple experiments greatly en- 
hanced our understanding of 613. In 201 1, through mea- 



surements of v e appearance from a beam, the long 
baseline experiments T2K [9] and MINOS [10] reported 
hints of a non-zero 613 at about 2.5 and 1.4 standard 
deviations 2 , respectively. In January 2012, the reactor 
antineutrino experiment Double-CHOOZ [11] also re- 
ported a hint of a non-zero 613 at 1.6 standard devia- 
tions with a single detector. The Daya Bay experiment, 
with six functionally identical detectors at three loca- 
tions, carried out a measurement of relative ratio of reac- 
tor antineutrino rates [12], which significantly improved 
the sensitivity to the sin 2 2613. In March 2012, the Daya 
Bay collaboration announced a non-zero value of 613 at 
5.2 standard deviations [13]. About one month later, this 
finding was confirmed by the RENO reactor antineutrino 
experiment [14], which reported a consistent result us- 
ing a ratio between antineutrino rates from two detec- 
tors. The Daya Bay experiment has since reported an im- 
proved measurement of the electron antineutrino disap- 
pearance [15] with 2.5 times of the previously reported 
statistics [13]. 

THE DAYA BAY EXPERIMENT 

The Daya Bay experiment, located on the south coast of 
China (55 km northeast to Hong Kong and 45 km east 
to Shenzhen), was designed to provide the most precise 
measurement of 613 with a sensitivity of sin 2 20i3 < 0.01 
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2 Results from the long baseline experiments actually depend on the 
assumption of the neutrino mass hierarchy and value of CP phase S. 
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FIGURE 1. The layout of the Daya Bay experiment. The 
dots represent reactor cores, labeled as Dl, D2, and Ll-4. 
Six ADs were installed in three experimental halls during the 
reported analysis period. 




FIGURE 2. Layout of Daya Bay detectors in a near site. 



at a 90% C.L. [16]. Such a measurement requires high 
accuracy and precision. The high accuracy is achieved 
by the combination of powerful reactors (17.6 GW ther- 
mal power) and large target mass (80 tons in the far 
hall). In addition, the location of the far detectors is op- 
timized to obtain the best sensitivity to sin 2 2 613 with 
the current knowledge of AmZ 2 - To achieve high preci- 
sion, the reactor-related systematic uncertainties are min- 
imized by adapting the ratio method [12] with multi- 
ple detectors at multiple baselines. The detector-related 
systematic uncertainties are minimized by using iden- 
tical detectors and performing precise detector calibra- 
tions. The background-related systematic uncertainties 
are minimized by placing detectors deep underground 
in order to reduce cosmic muon related backgrounds. 
Furthermore, passive shielding (water pools surrounding 
detectors) and active shielding (resistive-plate chambers 
above water pools) were implemented to tag the cosmic 
muons. The water pools also shield detectors from vari- 
ous environmental radioactive backgrounds. 

Fig. 1 shows the layout of the Daya Bay experiment 
during this reported analysis period. There are six re- 
actors grouped into three pairs. Each reactor contains a 
core with a maximum 2.9 GW thermal power. Three un- 
derground experimental halls (EHs) are connected with 
horizontal tunnels. The effective vertical overburdens 
are 250, 265, and 860 water-equivalent meters for EH1, 
EH2, and EH3, respectively. For this improved measure- 
ment, two, one, and three antineutrino detectors (ADs) 
were installed in EH1, EH2, and EH3, respectively. The 
distance from the six ADs to the six cores were surveyed 
with the Global Positioning System (GPS) above ground 
and Total Stations underground. The precision of dis- 
tance was about 1.8 cm. 

The antineutrinos in the Daya Bay experiment are 
detected through the inverse beta decay (IBD) process 



V e + p — >• e + + n. The signature of such process is a 
prompt signal of the scintillation and subsequent anni- 
hilation of the position in the liquid scintillator (LS), fol- 
lowed shortly by a delayed signal with —8 MeV energy 
deposition when the neutron is captured by the Gadolin- 
ium (Gd) doped inside LS (0.1% in weight). The energy 
of neutrino can be deduced from the energy deposition 



of the prompt signal with Ey ps E 



positron ' 



0.8 MeV. As 



shown in Fig. 2, the ADs adopt a three-zone cylindrical- 
shaped design, with the inner, middle, and outer zone 
containing 20 ton Gd-doped LS (Gd-LS), 20 ton LS, and 
40 ton mineral oil, respectively. With load cells and an 
ISO tank, the target mass uncertainty of the 20 ton GD- 
LS is controlled to be only about 3 kg. 

Each AD contains 192 8 -inch photomultiplier tubes 
(PMTs) installed on the side walls. The photocathode 
coverage is about 8%, which is further enhanced to about 
12% with a pair of optical reflectors at the top and bottom 
of each detector. The achieved detector energy resolution 
is parametrized as [17] 
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with respect to the visible energy E. The detector calibra- 
tion is performed weekly with three automated calibra- 
tion units (ACUs) per AD: two located above the center 
and the edge of the GD-LS region and one placed above 
the LS region. The ACUs are remotely controlled by a 
Lab VIEW program and operated underwater. Each ACU 
contains four sources: a light-emitting diode (LED) for 
the PMT gain/timing calibration, a — 15 Hz 68 Ge source 
for the threshold calibration of the IBD prompt signal, a 
-100 Hz 60 Co source for the high statistical determina- 
tion of the overall energy scale, and a —0.5 Hz 241 Am- 
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FIGURE 3. Discrimination of the flasher events (FID > 0) from the delayed signals of IBD (FID < 0). While the flasher events 
possess different distributions among different ADs, the delayed signals of IBD share the same distribution. 



C neutron source to understand neutron captures on Gd 
and to determine the H/Gd ratio in the target region. 

The muon detection system in each experimental hall 
consists of a high purity water pool surrounding the ADs 
and a layer of resistive plate chambers (RPC) above the 
water pool. The water pool is further divided into two 
optically separated regions, the inner water pool (IWS) 
and the outer water pool (OWS). Each region operates 
as an independent water Cerenkov detector, and can be 
used to cross calibrate each other. The muon detection 
efficiencies are measured to be 99.7% and 97% for the 
IWS and OWS [17], respectively. The water pool also 
plays a crucial role in shielding radioactive backgrounds. 
The distance between the edge of each AD to the closest 
wall is at least 2.5 m. The array of RPCs covering the 
entire water pool is used to provide additional tagging of 
cosmic muons. 



SELECTION OF INVERSE BETA DECAY 
EVENTS 

About 5% of the PMTs would spontaneously flash and 
emit light. Such events are called "flashers". The recon- 
structed energy of such events covers a wide range, from 
sub-MeV to 100 MeV. A few features are observed when 
a PMT flashes: i) the charge fraction of the flashing PMT 
is high; ii) the surrounding PMTs as well as the ones lo- 
cated on the opposite side of the AD receive large frac- 
tion of light from the flashing PMT; and iii) the timing 
spread of all PMTs' hits are generally wide. Accord- 
ingly, a few flasher identification (FID) variables were 



constructed to separate the good physics events from the 
flashing events. Fig. 3 shows the distribution of an FID 
variable deduced from the charge pattern for the IBD 
delay candidates. The good IBD events are well sepa- 
rated from the flasher events. Detailed description of the 
flasher discrimination can be found in Ref. [17]. The in- 
efficiency of the IBD signals due to FID cuts is only 
about (0.024 ± 0.006)%, and the contamination of flash- 
ing events in the IBD sample is below 0.01%, which is 
further suppressed by the accidental background subtrac- 
tion procedure. 

After the flashing events are removed, the IBDs are 
further selected with the following cuts: i) the energy of 
the prompt signal is between 0.7 and 12 MeV; ii) the en- 
ergy of the delay signal is between 6 and 12 MeV; and 
iii) the time difference between the prompt and the delay 
signal is between 1 and 200 jJ.s. In addition, a multiplic- 
ity cut is applied to remove the energy ambiguities in the 
prompt signal. For example, one choice of the multiplic- 
ity cut requires no prompt-like signal 400 /i s before the 
delay signal and no delay-like signal 200 jj.s after the de- 
lay signal 3 . The fixed time cut (relative to the delay sig- 
nal) leads to simplified calculations of the efficiency and 
livetime of IBD events and the rate of accidental back- 
grounds. Correspondingly, three types of muon veto cuts 
are also applied to the delay signal in order to suppress 
backgrounds. The first one is for the water pool muon, 
which is defined as one IWS or OWS event with more 



3 The prompt-like and delay-like signals refer to events with energy 
0.7-12 MeV and 6-12 MeV, respectively. 



than 12 PMT hits. The veto cut spans from 2 /is before to 
600 /is after the water pool signal. The second one is for 
the AD shower muon, which is defined as one AD signal 
with more than 3 x 10 s photoelectrons (PEs). The corre- 
sponding energy is about 1.8 GeV. The veto cut spans 
from 2 /is before to 0.4 s after the AD signal. The last 
one is for the AD non-shower muon, which is defined as 
one AD signal with its energy between 20 MeV and 1.8 
GeV. The veto cut spans from 2 /is before to 1.4 ms after 
the AD signal. The choice of 2 /is before the muon signal 
in all three cuts is to leave enough room for the potential 
non-synchronization among different detectors. The AD 
shower muon veto is applied to suppress the long-lived 
9 Li/ 8 He background. The AD non-shower muon veto is 
applied to suppress the fast neutron background. The 
overall detection efficiency for the IBD events is about 
80%, with the Gd catpure ratio (84%) and the efficiency 
of the 6 MeV delay signal cut (91%) being the two lead- 
ing contributors. 

SUMMARY OF BACKGROUNDS 

The largest contamination in the Daya Bay IBD sam- 
ple is the accidental background with about 4.6% at the 
far site and about 1.7% at the near sites. An accidental 
background event arises when a delay-like signal and an 
unrelated prompt-like signal (usually radioactive events) 
accidentally fall within a 199 /is coincidence window. 
Such background {R aC cidentai) can be accurately calcu- 
lated through Poisson statistics given the measured rates 
of single prompt-like events (R p ) and single delay-like 
events (R c i) 4 : 

Raccidental = P(0, 200/U • R p ) • P(l , 199/1.9 • R p ) 

■R d P(0,W0flsR d ). (2) 

Poisson function P(n,fl) = represents the prob- 

ability of observing n events given an expectation value 
of /i events. The above calculation results in negligible 
systematic uncertainties. This method is cross-checked 
with an off-window coincidence method and a coinci- 
dence vertex method [15]. 

The second largest contamination at the far site is the 
correlated backgrounds from the Am-C neutron source. 
During the data taking, the Am-C neutron sources are 
parked inside the ACUs on top of the ADs. The energetic 
neutrons from these sources occasionally undergo an in- 
elastic scattering with an iron nuclei resulting in gamma 
emissions, followed by a neutron capture on another iron 
nuclei with additional gamma emissions. When these 



4 The efficiency due to the fixed-time multiplicity cut is then 
P(0, A00\ls • R p ) • P(0, 200ns ■ R d ) . 




FIGURE 4. Three ADs were deployed in the EH3. 



gammas are emitted toward AD, it is possible these cor- 
related events could mimic an IBD signal. The contam- 
ination from the Am-C source, which is about 0.3% 
(0.03%) at the far (near) site, is calculated through a 
GEANT4-based Monte Carlo simulation, which can re- 
produce the energy spectrum of single backgrounds from 
Am-C sources. The systematic uncertainty of this corre- 
lated background is assumed to be 100%. 

The 9 Li/ 8 He and the fast neutrons are two major IBD 
contamination caused by cosmic muons. The nucleus of 
9 Li or 8 He are produced from the carbon nucleus when 
cosmic muons pass through the liquid scintillator. The 
9 Li and 8 He, with half lifetime of 257 ms and 172 ms, 
respectively, are both long lifetime beta emitters. They 
would undergo a beta decay providing a prompt-like sig- 
nal. The daughter nuclei could then undergo a sponta- 
neous fission, with a neutron emission in the final state 
resulting in a delay-like signal. Such a pair of corre- 
lated prompt-like and delay-like signal would mimic an 
IBD event. The contamination of 9 Li/ 8 He can be directly 
measured by fitting the spectrum of time between the 
IBD candidate and the last tagged AD muon. The mea- 
sured 9 Li/ 8 He rates from all three experimental halls are 
consistent with an empirical formula of OCE^ 74 given 
the average muon energy in each site. Furthermore, 
the contamination are suppressed by an optimized AD 
shower muon cut described in the previous section. The 
remaining contamination is about 0.2% (0.35%) for the 
far (near) site with a 50% systematic uncertainty due to 
the fitting procedure. 

The fast neutron backgrounds are caused by the ener- 
getic neutrons produced inside or outside the muon veto 
system. These energetic neutrons can undergo an elas- 
tic scattering with protons, leaving a prompt-like signal 
due to the proton recoil, followed by the neutron ther- 
malization and then neutron capture on Gd producing a 
delay-like signal. The energy of the proton recoil signals 
ranges from sub MeV to tens of MeV. Therefore, one 
can extrapolate the measured fast neutron's prompt en- 
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FIGURE 5. The ratio between the measured and the expected signals in each detector vs. the flux-weighted average baseline, 
which is computed with reactor and survey data. The expected signals have been corrected with the best fit normalization parameter 
assuming no oscillation. The error bars represent the uncorrected uncertainties. The oscillation survival probability at the best-fit 
value is plotted as the smooth curve. For clarity purpose, the AD4/6 data points are shifted by -30 and 30 m, respectively. The y} 
value vs. the sin 2 20i3 value is shown in the inner panel. 



ergy spectrum above 15 MeV to the energy range of in- 
terests (0.7-12 MeV) in order to estimate the contamina- 
tion in the IBD candidates. A fiat background spectrum 
is assumed, which is confirmed by the spectrum of fast 
neutron events with muon tagging from water pools and 
RPC. The systematic uncertainty is constrained to about 
30%. The contaminations due to fast neutron are esti- 
mated to be about 0.07% (0.12%) for the far (near) site. 

The last contamination is the 13 C(a,n) 16 background 
caused by radioactivity inside ADs. The contamination is 
determined by Monte Carlo with measured alpha-decay 
rates. We identified four sources of alpha decays: the 
210 Po events and the decay chains from 238 U, 232 Th, and 
227 Ac. The backgrounds are calculated to be about 0.05% 
(0.01%) for the far (near) site with a 50% systematic 
uncertainty. Altogether, the total backgrounds in the IBD 
sample are thus determined to be 5 ± 0.4% and 2 ± 0.2% 
for the far site and near sites, respectively. 

OSCILLATION ANALYSIS 

The oscillation analysis includes data of 139 days (Dec. 
24th 201 1 - May 1 1th, 2012) with six ADs. Fig. 4 shows 
a picture of the three ADs in the EH3. The analysis pro- 
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FIGURE 6. Top panel: Measured prompt energy spectrum at 
the far site (sum of three ADs) is compared with the predicted 
spectrum based on the measurements at the two near sites 
assuming no oscillation. Backgrounds are subtracted. Only 
statistical uncertainties are shown. Bottom panel: The ratio of 
the measured over the expected (no-oscillation) spectrum. The 
solid curve is the expected ratio vs. the true prompt energy with 
the neutrino oscillation for sin 2 20i3 = 0.089. 



cedure is described in details in Ref. [13, 15]. Since dif- 
ferent experimental halls have different mountain over- 
burdens, they also observe different muon rates. The 
average muon veto efficiency for six ADs in the three 
experimental halls are 0.8231, 0.8198 (EH1), 0.8576 
(EH2), and 0.9813, 0.9813, 0.9810 (EH3). Besides the 
background-related uncertainties, the largest uncorre- 
lated detector-related uncertainty (0.12%) is due to the 
6 MeV energy cut in selecting delay signals. Other siz- 
able uncorrelated detector-related uncertainties include 
0.1% from the neutron Gd capture ratio, 0.03% from the 
number of target protons, and 0.02% from the spill-in 
effect. (The spill-in effect refers to that the IBD neu- 
trons generated outside but captured inside the target 
GD-LS region outnumber the IBD neutrons generated 
inside but captured outside the GD-LS region. The rea- 
son for such imbalance is that thermal neutrons have a 
larger cross section to be captured on the Gd than the 
proton.) With the ratio method, the correlated detector- 
related uncertainty (about 1.9% in total) has negligible 
effects on the oscillation analysis. The same applies to 
the correlated reactor-related uncertainty. The uncorre- 
lated reactor-related uncertainties include 0.5% from the 
received thermal power data, 0.6% from the calculated 
fission fractions, and 0.3% from the antineutrinos pro- 
duced by the spent fuel. The total uncorrelated reactor- 
related uncertainty is 0.8%, which is further suppressed 
by about a factor of 20 in the oscillation analysis due to 
the multiple core/reactors configuration at Daya Bay. 

The final antineutrino rates per day in the six ADs 
of three experimental halls, after corrections of the live- 
time, the veto efficiencies, and the background, are 
662.47±3.00, 670.87±3.01 (EH1), 613.53±2.69 (EH2), 
and 77.57±0.85, 76.62±0.85, 74.97±0.84 (EH3). In the 
same experimental hall, the differences of AD rates stem- 
ming from ADs' different locations are consistent with 
expectation within statistical uncertainties. The deficit of 
antineutrino rate at the far site is quantified by the ra- 
tio of the measured far-hall IBD rate over the expected 
rate, which is calculated with the measured IBD rates of 
the near detectors assuming no oscillation. The resulting 
deficit is 0.056 ± 0.007 (stat.) ± 0.003 (sys.). 

The chi-square method with pull terms is used to 
extract sin 2 20i3 within the standard 3-flavor oscillation 
model, in which the disappearance probability of the 
electron antineutrino is written as: 

P(v e — > V e ) = 1 —sin 2 2 0i 3 cos 2 0i2 sin 2 

- sin 2 20i3 sin 2 0i2 sin 2 A 32 

— cos 0n sin 20i2 sin A21, (3) 
u, 
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the uncertainties from the backgrounds, detectors, reac- 
tor fluxes, and the oscillation parameters are taken into 
account properly in a consistent manner. The sin 2 20i3 
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FIGURE 7. Projected sin 2 20i3 's sensitivity in terms of half 
of the 68% confidence interval of Daya Bay vs. the running 
time. The published and improved Daya Bay results are shown 
as red dots. The period of 6- AD data taking (black) is sepa- 
rated from the period of 8- AD data taking (blue) by an installa- 
tion/calibration period (dashed). 



is determined to be 0.089±0.011. Fig. 5 shows the ra- 
tios of measured over expected IBD rates vs. weighted 
baseline for all ADs. The data are compared with the ex- 
pected oscillation curve (red-solid line). Our improved 
measurement disfavors the sin 2 20i3 = hypothesis at a 
7.7 standard deviations. 



CONCLUSION AND OUTLOOK 

In this talk, we reported an improved measurement of 
electron antineutrino disappearance at Daya Bay. The 
new results with 2.5 times more data confirms the pre- 
viously published results [13], and improves the preci- 
sion of sin 2 20i3. Currently, the total uncertainty is still 
dominated by the statistical uncertainty, which is about a 
factor of 2 larger than the systematic uncertainty. Fig. 7 
shows the projected sensitivity of sin 2 20i3 in terms of 
half of the 68% confidence interval with respect to the 
running time. With the full 8 -AD configuration, we ex- 
pect to achieve a ~ 5 % measurement of sin 2 2 0i 3 in about 
3 years. Furthermore, the Daya Bay experiment also has 
the potential to measure the effective squared-mass dif- 
ference Am 2 e , which is a combination of Am 2 j and Am 2 2 , 
through the measurement of IBD energy spectrum distor- 
tion. Due to the short baseline < 2 km, the measurement 
of Am 2 ,, is not sensitive to the neutrino mass hierarchy 
(sign of Am 2 2 ) [18]. In addition, the high statistics IBD 
samples from Daya Bay would provide the most precise 
measurement of the antineutrino energy spectrum, which 
is essential for the future measurements with reactor an- 
tineutrinos. 



With the global efforts led by the Daya Bay exper- 
iment, the sin 2 20i3 is found to be around 0.09. Such 
a large value of sin 2 2 613 opens doors to two of the 
remaining unknown parameters in the neutrino sector, 
the CP phase 8 in the leptonic sector and the neutrino 
Mass Hierarchy 5 . In particular, the long baseline exper- 
iments [20, 21, 22, 23] can provide essential information 
for both parameters through the (anti) v e appearance from 
a (anti)v^ beam. Meanwhile, the possibility of utilizing a 
medium baseline (^60 km) reactor neutrino experiment 
to determine the neutrino mass hierarchy is also inten- 
sively discussed [24, 25, 26, 27, 18, 28]. We therefore ex- 
pect a new era of discovery in the next couple of decades. 
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A Bayesian approach in presenting results/sensitivities of the neu- 
trino mass hierarchy has been proposed in Ref. [19]. 



